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Abstract 


This report describes an application of the method of moments to 
calculate the electromagnetic scattering from irregularly shaped , thin, 
metallic flat plates in free space . In the present technique , an ir- 
regularly shaped plate is enclosed by a rectangle on which the surface- 
current density is then expressed in terms of subdomain functions by 
dividing the rectangle into subsections. A shape function is introduced 
to ensure zero current outside the patch. The surface- current density is 
determined using the electric field integral equation (EFIE) approach 
in conjunction with the method of moments, and from a knowledge 
of the surface- current density , the electromagnetic scattering from a 
plate is calculated. Using this technique, the electromagnetic scattering 
from (1) a hexagonal plate, (2) an equilateral triangular plate, (3) an 
equilateral triangular plate with a concentric, equilateral triangular hole 
and an inverted , equilateral triangular hole, and (4) a diamond- shaped 
plate is computed and compared with the numerical results obtained by 
using the Electromagnetic Surface Patch (ESP) code developed by Ohio 
State University. The numerical results compare favorably with the 
measurements performed on these shapes in the Langley Experimental 
Test Range facility . 


Introduction 

A knowledge of electromagnetic (EM) scattering 
from a complex-shaped metallic object is of practical 
interest to electromagnetic analysts and engineers. 
In recent years, considerable interest has developed 
in the EM scattering analysis of polygonal plates 
because a complex-shaped metallic object can al- 
ways be modeled as an interconnection of these poly- 
gonal flat plates (ref. 1). An EM scattering anal- 
ysis of a polygonal, metallic flat plate can be 
accomplished by using various numerical techniques. 
One of the widely used techniques is the method 
of moments in which a polygonal plate is first seg- 
mented into a number of rectangular or nonrectangu- 
lar surface patches called subdomains. The unknown 
surface-current density on these subdomains is then 
determined using the electric field integral equation 
(EFIE) in conjunction with the method of moments. 

The early work on polygonal flat plates (ref. 2) 
consisted of subdividing an irregularly shaped plate 
into a collection of curvilinear cells that were rep- 
resented numerically as fourth-order polynomials. 
However, use of that method required a knowledge 
of the coordinates of a large number of points on the 
object, and this made the method unattractive. In 
another approach (refs. 3 and 4), a nonrectangular 
plate was viewed as an interconnection of polygonal 
plates. The unknown amplitudes of the surface patch 
modes were determined by expressing the currents on 
the polygonal plates in terms of nonrectangular sur- 


face patch modes and using the EFIE in conjunction 
with the method of moments. 

An Electromagnetic Surface Patch (ESP) code 
was developed in 1988 under NASA Grant NSG 1498 
by E. H. Newman at the ElectroScience Laboratory, 
Ohio State University. (This code is described in 
a manual entitled A User 7 s Manual for the Electro- 
magnetic Surface Patch Code: ESP Version IV.) The 
ESP code was based on the segmentation techniques 
described in references 3 and 4. Even though the 
segmentation technique used in the ESP code com- 
pletely filled the area of the nonrectangular plate, 
extra current modes were required at the intercon- 
nections of the polygonal plates to ensure a conti- 
nuity of current. Furthermore, when nonrectangular 
surface patch modes were used to express the sur- 
face current, the resulting matrix would be symmet- 
ric but not of a block Toeplitz nature, and hence it 
would require more impedance matrix filling .time. 
In the present work, a segmentation technique that 
is both symmetric and results in a block Toeplitz 
impedance matrix is developed to analyze the EM 
scattering from polygonal, thin, metallic flat plates. 

In the present analysis, an irregularly shaped 
plate is assumed to be enclosed by a rectangle with 
sides equal to W x and W y , the maximum dimen- 
sions in the x - and y-directions, respectively. By 
dividing W x into (M + 1) subdivisions and W y into 
(N 4-1) subdivisions, the surface-current density over 


the rectangle is expressed in terms of overlapping 
triangular functions in the direction of current flow 
and a pulse function in the orthogonal direction. Zero 
current outside the plate is ensured by introducing a 
space function in the current expansion function. (A 
space function is equal to 1/0 if the subdomain lies 
inside/outside of the irregular plate.) Selecting the 
testing functions to be the same as the expansion 
functions allows the EFIE to be reduced to a matrix 
equation that is solved by using standard matrix- 
equation solver subroutines. The surface. current on 
the plate is then used to determine EM scattering due 
to the plate. By using this technique, the EM scat- 
tering due to several nonrectangular plates is com- 
puted and compared with the results obtained from 
Newman’s ESP code. The results obtained from us- 
ing the present method are also compared with the 
experimental data measured in the Langley Experi- 
mental Test Range (ETR) facility, which is a compact 
range specifically designed for microwave scattering 
measurements. The ETR is a dual anechoic chamber 
with a Greorian reflector system (ref. 5) containing 
a 16-ft 2 rolled-edge main reflector. 

The measured data are provided on a floppy disk 
for the reader, a description of which is given in the 
appendix. 

Symbols 


A (x,y, z) 

magnetic vector potential 

E * 

incident electric field vector 

|Ej| 

absolute value of incident electric 
field vector 

E s [x,y,z) 

scattered electric field vector 

E SX 

x-component of scattered field E s 

E S y 

y-component of scattered field E s 

E s0 

^-component of scattered field E s 

E S<fi 

(^-component of scattered field E 6 - 

E xi,Eyi t E z i 

x-, y- 3 and ^-components, respec- 
tively, of incident electric field 

E <p, 

0- and (^-components, respec- 
tively, of incident electric field 

e i(^X! Ay , 0) 

Fourier transform of E sx (x,y, 0) 

e y (k x , k y , 0) 

Fourier transform of E sy (x, y,0) 

E xmn(kx i ky) 

Fourier transform of 
Pjn(x) Qn{y ) 


E ymn{kxi ky) 

Fourier transform of 
Qm{ x ) Pn{y) 

{FX} 

sin (k x Ax/2) 
~ k x Ax/2 

{FY} 

_ sin(fc v Ay/2) 
ky Ay /2 

{FX1} 

_ sin[(sin 0 t cos <^ z )(^0 Ax/2)] 
~ sin 9i cos <j>i ko(Ax/2) 

{FYl} 

_ sin [(sin $i sin M(ko Ay/2)} 
(sin 0 l sin (j>i)(k{) Ay/2) 

/ 

frequency, GHz 

G(x,y, z/x' ,y') 

free-space- scalar Green’s function 

g(k x , k y , z/x', y') Fourier transform of G{x, y, zjx f , y 1 


incident magnetic field vector 

h(p) 

complex amplitude of pth x- 
directed subdomain current mode 

hio) 

complex amplitude of qth y- 
directed subdomain current mode 

J s(x.y) 

induced surface-current density 
vector on plate 

Jsx{% j y) 

x-component of J 5 (x,y) 

Jsy(x, y) 

y-component of J 6 (x,y) 

j 


j.s 

Fourier transform of J s 

jsx 

x-component of } s 

Jsy 

y - component of j iS 

k* 

propagation vector of plane wave 

kx 

Fourier transform variable with 
respect to x 

ky 

Fourier transform variable with 
respect to y 

k z 

complex propagation constant in 
z-direction 

ko 

propagation constant in free 
space 

M -j- 1 

number of subdivisions in 
x-direction 

m, n 

(m, n)th subdomain of induced 
current 

m 7 , n f 

(ra 7 , n 7 )th subdomain of induced 
current 

N+ 1 

number of subdivisions in 
y-direction 
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p 

total number of x-directed sub- 
domains on plate 

Pm{ x ) 

piecewise linear distribution in 
x-direction 

Pn{y) 

piecewise linear distribution in 
y-direction 


/ 

Zxy mutual impedance between qth 

y-directed and p f th x-directed 
subdomain currents 

Zyx mutual impedance between pth 

x-directed and q f th y-directed 
subdomain currents 


P 

equivalent to (m, n)th x-directed 
subdomain 

* yy 

P f 

equivalent to (m f , n^th x- 
directed subdomain 


Q 

total number of y-directed sub- 
domains on plate 


Qm{y) 

pulse distribution in y-direction 

ot 0 

Qtl (*^) 

pulse distribution in x-direction 

Ax 

Q 

equivalent to (m, n)th y-directed 
subdomain 

A, 


mutual impedance between qth 
and q'th y-directed subdomain 
currents 


polar coordinate variables related 
tO kg } ky 

angle between Ej and 6 deg 


W x 

7f7 + T 


Wy 

W+1 



equivalent to (m / ,n / )th y-directed 
subdomain 

R p J 

real part of Z X x 

V 

position vector in direction of 
plane wave 

S x (m,n) 

shape function for x-directed 
current (= 1 or 0) 

S y (m , n ) 

shape function for y-directed 
current (= 1 or 0) 

v/ 

reaction of p'th x-directed sub- 
domain testing function with E x i 

y/ 

reaction of q* th y-directed sub- 
domain testing function with E y i 

W x 

maximum dimension of plate in 
x-direction 

Wy 

maximum dimension of plate in 
y-direction 

yP'P' 

imaginary part of Z x f 

z,y,z 

Cartesian coordinates of field 
point 

x\y\z f 

Cartesian coordinates of source 
point 

^xx 

mutual impedance between pth 
and p' th x-directed subdomain 
currents 


6 delta function 

770 free-space wave impedance 

cf> angles of scattered electro- 

magnetic wave 

9 i,<f>i incident angles of electromagnetic 

wave 

0 n fa unit vectors along 61 and (j>{ 

axes, respectively, in spherical 
coordinate system 

0 S , (j) s unit vectors along 6 S and <j) s 

axes, respectively, in spherical 
coordinate system 

Ao wavelength in free space 

po magnetic permeability of free 

space 

a total radar cross section 

<7 he copolarized radar cross sec- 

tion when E-polarized wave is 
transmitted 

c 7 eh cross-polarized radar cross 

section when E - polarized wave 
is transmitted 

<7 HE cross-polarized radar cross 

section when H - polarized wave 
is transmitted 
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&HH 

copolarized radar cross sec- 

Abbreviations: 



tion when H - polarized wave is 
transmitted 

EFIE 

electric field integral equation 



EM 

electromagnetic 

LJ 

angular frequency, 2nf 

ESP 

Electromagnetic Surface Patch 

V 

gradient operator 

RCS 

radar cross section 

Theory 





General Theory 

Consider an irregularly shaped, infinitesimally thin plate excited by a plane wave as shown in figure 1(a), 
By using the representation of reference 6, the incident field with a time variation of e Jujt may be written as 

Ei(ar, y, z) = + frE#) e~i k > T 

= (&i |E,;| cosQo+^i |Ej| sin a 0 ) e“^ k, r (1) 


where 


kj • r = — fco sin 9 j (x cos (f>i + y sin <f> j) 

and Uq is the free-space wave number. From equation (1), the x-, y-, and 2 -components of the incident field 
may be written, respectively, as 

E x i - E 0 . cos 9 t cos 4>i - Eq. sin 4> l 

E yi ~ E 0, cos & i sin & + e 4>, cos fc ( 2 ) 

E z i = ~Eg i sin 9 8 


The corresponding magnetic field components are obtained through 


Hi = — k t x E i 

% 


( 3 ) 


where 7 ^ is the free-space impedance. The incident field with E$ % ^ 0, E^. = 0 (i.e., ao = 0°) is called the 
H-polarized wave and with Eq. ~ 0, E^. ^ 0 (i.e., oq = 90°) is called the E-polarized wave . 

Let Js(:r,y) be the induced surface-current density on the plate. The electromagnetic field due to J iS (:E,y) 
located in the z = 0 plane may be obtained from the magnetic vector potential A (x,y,z) as 


H* (x, y, z) = — - V x A(x, y, z ) (4) 

Mo 

E s (x, y, z) = |/cq A(x, y, z) + V [V ■ A(x, y, 2 )]| (5) 

where //q is the permeability of the medium. The vector potential A(x, y, z) in equations (4) and (5) satisfies 
the wave equation 

V 2 A(x, y, z) + hi A(x, y, z) = -p 0 y) (6) 

If G(x y y, zjx\ y f ) is the free-space-scalar Green’s function, the magnetic vector potential A(x, y, z) is obtained 
from 

A(x, y, z) = J J J s (x',y') G{x,y,z/x' ,y')dx' dy ' (7) 

where G(:r, y, zjx ! , y r ) satisfies 


(v 2 + A: 2 ) G(x, y, zjx\ y) = -^o ~ *’) 6{y - y) 6{z) (8) 
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The solution of equation (8) in the (k x , k y ) domain may be written as (ref. 7) 


OO OC 

G(x,y,zlx',y')=^ J dk, J g(k„ k„, z/x', y' ’) dk 


( 9 ) 


k x —~ OC ky=— 


where 


g(.kx,ky,z/x',y') - 


k, = 


y/ko- k l~ k v ( k l + k l^ k o) 

-j \f k x + - kl (kl +k 2 y > kl) 


and the + and - signs in the exponential are used for z < 0 and 2 > 0, respectively. Substituting equation (9) 
into equation (7) gives the magnetic vector potential 


A(x ’ 2/ ’° ) = 


OO OO 

J dkx J e jk x x+jk y y 


dk,, 


k x =—oc ky —— oo 


where 


}s( k x, ky) = xjsx + yjsy 

= f dx' J J s(x',y') e~i kxX '~i kyy ' dy' 


( 10 ) 


( 11 ) 


Substituting equation (12) into equation (5) allows us to write the scattered tangential components of the 
electric field over the plate due to the induced currents as 


DO OO 

E S x(x,y, 0 ) = (27t/ 2 P I ( ^ x I [(^0 _ k x) jsx(k x ,ky) — k x ky j 3 y(k x ,ky)^ xe? kx 
^ k T ——oc k v =— oo 


x+jk y y dky ( 19-) 

2 k z K ’ 


E S y(x, y, 0) - 


OC OO 

J dk x J ^( — k x ky) jsx(k x , ky) + ~ ky^j jsy(k x ,t 


k x =— oo ky— — oo 


(13) 


x e jkxX+jkyy dky 

2k z 


Subjecting the total tangential electric field on the plate to zero gives the following equations 

E sx + E xt e^ k ‘ r) = 0 
E sy + E yi e-^ T) = 0 

Substituting equations (2), (12), and (13) into the above equations gives the following integral equations with 
surface current as an unknown variable: 


(2tt) 2 ^q 


OO oo 

j dk, J [(«$-*») i„(k„k,)-k x ky i„(k„ *„)]* e> k ‘^y = E xi e-i*' (14) 


k x =— oo k y =—rx> 
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(2 


OC DC 

^2 J dk * J Js,(k x , ky) + (fc 0 2 - k 2 y ) j sy {k x , *„)] x e ikxX+jk y y ^ = £^. e -JkiT (15) 


k x ~—o O ky=—0 o 


With the aid of figure 1(b), the surface- current density distributions can be expressed as 

M N + 1 

Jsx(x,y) - E 53 s x( m > n ) Ix{m,n) P m {x) Q n {y) 
m~ 1 n=l 


and 


M+l A r 

^sy( x > 2/) = EE Iy(m,n ) Q m (x) P„(j/) 

m=l n— 1 


where /^(m, n) and /^(m, n) are the unknown current amplitudes and 

P m{ x ) ~ | 


1 

to 


Qn{y) — 


f i" *£ 

x— X| 


x 


Ax 


((arm - Ax) < x < x m) 
fam < x < (x m + Ax)) 
(Otherwise) 


(n — 1)A y < y < n Ay) 
(Otherwise) 


If the (m,n)th cell lies inside the plate, then 


S x (m,n) 1 = 
S y {m,n) j 


1 


Otherwise, 


Also, 


S x {rn, n) \ __ 
5„(m, n) J “ 


Ax = 


Wr 


M + 1 
= ^ 


(16a) 


(16b) 


N - hi 

Here, and W y are the maximum dimensions of the plate in the x- and y-directions, respectively. The 
(m, n)th cell is considered to be inside the plate if the area occupied by the plate is more than 50 percent of 
the cell area. 

Using equation (11) allows us to write the Fourier transform of the patch current as 

M AT+1 

jsx{k x ,ky ) = E E 5 *( m > n ) * 

m—l n= 1 
M+l AT 

jsy(k x ,ky) = E E S y( m > n ) Jyirn^n) F 

m= 1 n=l 


l(^x » ky) 

(17a) 

(k # 5 ky) 

(17b) 
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where 


Fxmnikx, ky) = Ax A y {FY} {FX} 2 exp jk x (^y- - x m ) + jk y (^y - y n + ^y 

9 r fW x Ax\ (W v M 

Pymn(k Xl ky) — Ax {-^^0 ex P jkx ^”2 ^ 2~ J \/~ 2 l J 

sin ( k x Ax/2) 


{FX} = 
{FY} 


k x Ax/2 
sin (ky Ay/2 ) 


ky Ay/2 

With reference to figure 1(b), the double summation with respect to m and n in equation (17) can be 
represented by a single summation with respect to p or q. If P and Q are the maximum numbers of x- and 
y-domain subcells, respectively, on the plate, then equations (17) may be rewritten as 


jsx(kx) ky) — ^ ^ I x (p) F X p(kxi k y ) 

p= 1 

Q 

jsy(k x , ky) = lyio) Pyq{k X ) ky) 

<7=1 

Substitution of equations (18) into equations (14) and (15) and use of the method of moments yields 

P Q 

E 7 -(p ) z ™' + E z'F = v£ 

p=i 9=1 


(18a) 


(18b) 


(19) 


Q 


E 7 -(P) ZV J + E l M z ty = V v 

p = 1 9=1 

where p 1 = 1,2, ... P, q' = 1, 2, . . . Q, and 


( 20 ) 


yp,p f 

** XT ~ 


J J ^ 2 k - Fx ^ky) F xp /(k x , ky) dk x dky 


(27 r) 2 fcn , 


( 21 ) 


OO OG 

zlf = 7^2 / / Kp(k x ,ky) dk x dky (22) 

^ = — oo ky = — oo 

OO oo 

= ^27r)2fc2 J J Fxp(k x , k y ) F* q) {k x , k y ) dk x dk y ( 23 ) 

^ k x =-oc k y =-o o 

oo oo _ ^,2 

— ^^2 J J 2k F yq f(k x ,ky) dk x dk y ( 24 ) 

^ k x =—o o ky — oo 

V/ = J J E xl e~ jk ‘ T P m ,(x') Q n ’(y ) dx' dy ( 25 ) 

V/ = J j E yi e~ jKT Q m ,{x) P n i(y) dx' dy' ( 26 ) 
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In equations (21) (24), the superscript * indicates a complex conjugate. Using the expressions for F xp and F yq 
with the substitutions k x = ko0 cos a and k y - k Q 0 sin a allows equations (21), (22), (23), and (24) to be 
written, respectively, in terms of the new variables (a, 0) as 


y _ r A 2 x A 2 y 


Z p ’ p = 

" XX 


7T* 


oc n/2 

/W 


1 — 0 2 cos 2 a 


{FY} 2 {FX) A 


0=0 a=0 

x cos[fco(Xpp/) 0 cos a] cos^oft/pp/) 0 sin a] da 


where x pp t = — x m / and 


Vpfa ~ Vn Un') 


(27) 


z q,p 

^xy 


% k o a 2 z A2 y 


oc tt /2 

/ W- 


0 2 cos q sin a 

2vT^ 


{FU} 3 {FX} 3 


0=0 ct— 0 

x sin[fco(x ? y) (3 cos a] sin[fco(y 9 y) 0 sin a] da 


(28) 


where x^y x m x m! 3~ A j and y^y — ,Vn J/ ?l / 2 ^ > 


Z p ’ 9 ' = z Q ’ pl 

** yx ^xy 


and 


yy 


00 ^r/2 

/W 


y _ T7 qfcg A 2 x A 2 y f a Ja f 1 - 0 2 sin 2 a 

mt n 

7T Z 

0=0 a=0 


2\A - W 

x cos[fco(x 9 g/) 0 cos a] cos[fco(y g y) 0 sin a] da 


{ FYp{FX } 2 


(29) 


(30) 


Vq,q' — Vn - y n ' ■ In deriving equations (27) (30), the even and odd properties of 


where x^y = x m — x m t and 
the integrands have been utilized. 

Integrating equations (25) and (26) with respect to x’ and y' gives, respectively, 


Vx = Ax Ay E xi {FK1} {FX1} 2 exp 


jh 0 X 




sin &i cos fa 


x exp 


jk 0 y. 


Wy Ay 
~2 2 " 


sin 9i sin fa 


(31) 


and 


V q = Ax Ay E yi {FY l} 2 {FX1} exp 


.. ( W x Ax\ . 

JH) U m / — I sm Vi cos fa 


x exp 


i*o ( y n f — y ) sin ° l sin ^ 


(32) 


where {FX1} is obtained from {FX} by replacing 0 cos a with sin 6 j cos fa, and {FY 1} is obtained from 
{FY} by replacing 0 sin a with sin 6 l sin fa. 

The elements Z , Z%y , etc., of the coefficient matrix are determined from equations (27)-(30) by using 
gauss-quadrature numerical integration techniques. The presence of y / 1 - 0 2 in the denominators of the 
integrands of equations (27) (30) causes numerical difficulty around 0 = 1. However, this is avoided by 
splitting the 0 integration into two parts: the first part consists of integration from 0 = 0 to 1, and the second 
part consists of integration from 0 = 1 to oo. When 0 is replaced by sin £ in the first part and by cosh £ in 
the second part, equation (27) becomes 

Z p J = R™'+jX p J (33) 
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where 


„ 1.2 a 2„ a 2. V 2 V 2 

fig/ = ^ ° 2(7r)2 / smCdC J ( 1- sin 2 Ceos 2 a) {FY} 2 {FX} 4 

C-0 c*=0 

xcos[fco(^p,p') sin C cos a] cos[fco(T/py) sin £ sin a] da 
in k 2 A 2 r A 2 ?/ ^ 

*?/ = ° 2(7r)2 J cosh C rfC J (1 - cosh 2 C cos 2 a){FF} 2 {FX} 4 

C=0 a=0 

xcos[fco(xpy) cosh £ cos a] cos[fco(?/py) cosh £ sin a] da 

where { FX } and {FT} in the equations for R x f an d ^xx are J as defined earlier, with the appropriate change 
of variables. 

Similarly, the singularities due to the presence of \J 1 — /3 2 in the denominators of equations (28) -(30) are 
removed and thus yield expressions analogous to equation (33). 

Equations (19) and (20) can now be solved for I x (p) and I y {q) either by using the matrix inversion method 
or by decomposing the impedance matrix into lower and upper triangular matrices. We show later that the 
method of triangular decomposition (ref. 8) is preferred over the matrix inversion method because the former 
takes less central processing unit (CPU) time. 

Radar Cross Section 

The radiation field due to the plate is given by (ref. 9) 

Eds = [ e x{k x , k y , 0) cos 4> + e y (k x , k y , 0) sin 4>\ V (34) 

zirr 

E<ps = [ey{kx, k y ,0) cos 4> - e x (k x ,k y , 0) sin <j>] cos Q^-e~ : > kor (35) 

I'KT 

where e x (k x ,k y , 0) and e y (k x , k y: 0) are the Fourier transforms of the tangential electric field components at 
the interface z = 0 (obtained from eq. (13)) and are given, respectively, as 

—wllq [ (^0 — ^x) jsx(k x <, k y ) — k x k y j sy (k Xi k y )] 

ex(k x , k y , 0) = — p ± 1 (36) 

aJq 

and 

— ujllq | k x k y jsx{k x ,k y ) 4- (k^ — k y ) jsy{k x , k y )\ 

ey(k x , ky , 0) = —p. i — l (37) 

ruQ Z 

Substituting equations (36) and (37) into equations (34) and (35) gives, respectively, 

e S 9 = [?«*(**> M {!} + *»(**> M f 2 }] ^ e ~ jk ° r (38) 

and 

E s <fi = [jsx{kxi ky) {3} + jsy(kx , k y ) {4}] cos ^° r (39) 

where 

1,2 _ l 2 v, h, 

{1} = 2k z C ° s4> ~~2k7 Sm ^ 
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{2} = cos <P + 


kl - k 2 


2k, 


sin (p 


roi ~k x k y , 4 4-i 

{3} = ~2kf C ° S ^ _ 2 k7~ Sm ^ 


{4} = 


ju2 _ u2 

/VQ l\y 


kxk y 


cos cp + ■■•••- sin ^ 


2fc 2 ™ " ' 2k z 

By substituting k x = k ( j sin 0 cos <p and k y = fco sin 6 sin 0, the above equations reduce to 

{1} = — cos 4> cos 6 


{2} = sin cp cos 6 


{3} “ 
{4} = 


— /cq sin (p 
2 cos 6 
k() cos 4> 


2 cos 9 

By using the above expressions, the scattered electric field components are obtained as 
E s e = - 770 ^° S - [j sx (k x , k y ) cos <f> + j S y(k x , k y ) sin <f> ] ^e~ jk ° T 


(40) 


E s<p = -pr ky) cos 4> - jsx{k x , k y ) sin <p] jk[,r 


(41) 


After substituting for j x and j y and performing a few mathematical manipulations, the scattered field 
components are written as 


r, -j'nok o Ax fc 0 A y A 0 e i k ° T 
h s y = : := - cos w 


47t \fix 2 \fwr 


E T [ \ F rp{kx' ky) 

L A y COS * 

P= 1 y 


+ 


Q 

y 

9=1 


Fyqjkx ■ k y ) 
Ax Ay 


£/„(?) 


(42) 


^S(p 


-jr] 0 k 0 Ax fc 0 Ay A 0 e ^° T 


4^^777 2y/wr 


Q 


T,Vi) 


8=1 


ky) 

Ax Ay 


COS (p 


E T , \ Fxp(kx,ky) . 

Ax Ay 5m * 

P= 1 " 


(43) 


In equations (42) and (43), F xp (k x ,k y ) and F yq {k x ,k y ) are calculated at k x = -k 0 sin 6 cos 4> and 
ky = -ko sin 6 sin (p. The radar cross section of the plate, as defined in reference 10, is then obtained 
from 


where 


nlEJ 2 

a = lim 47rr — l0 
r_ *0° |E;| 2 


\E s \ 2 = \E se \ 2 + \E s<j> \ 2 
^ = \E 0i \ 2 + \E (j>i \ 2 


(44) 


The radar cross section defined in equation (44) is the total monostatic radar cross section (RCS) of an 
object. However, in most of the measurement, either ^-polarized or iJ-polarized waves are transmitted and the 
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/^-polarized and //-polarized scattered far fields are measured separately. In order to compare the calculated 
results with the measurement, the radar cross section may be defined as 


a HH = lim 4?rr 2 ^j 2 

r ^°° I e B i \ 2 

° HE = r-^oo ^ JE^I 2 
a EH = lim 47rr 2 j^ sg ] 2 

a EE = lim 47r r 2 

l^l 2 


(45a) 

(45b) 

(45c) 

(45d) 


Numerical Results 

In this section the RCS of nonrectangular, thin 
metallic plates is calculated by using equations (44) 
and (45). To validate the present formulation, the 
RCS computed by using equation (44) is compared 
with (1) the RCS computed by using Newman’s 
ESP code and (2) the RCS measured in the Langley 
Experimental Test Range facility. 

RCS of Hexagonal Plate 

As a first example, a hexagonal plate with a = 
2.074 cm, as shown in figure 2, is considered. To 
study the edge behavior, the induced current densi- 
ties along the XX- and YY-planes, as shown in fig- 
ure 2, arc calculated using the present method and 
are plotted in figure 3 for i?-polarized (uq = 90°) 
and H-polarized (»q = 0°) incident waves. Fig- 
ure 3 shows that the normal and tangential com- 
ponents of the current to an edge tend to behave 
as expected; i.e., the normal component of the cur- 
rent to an edge goes to 0 and the tangential compo- 
nent to an edge approaches oc. Figure 3 also shows 
that M = N > 19 gives stable values of the current 
densities. 

Furthermore, to establish the convergence of the 
method, monostatic RCS’s of a hexagonal plate using 
equation (44) are calculated as a function of angle 
of incidence 6 { for M ~ N = 15,19,21, and 23 
and are presented in figure 4. In figure 4(a), 0° < 
6 t < 90° for fa = 0° or 180° indicates the XZ- 
plane. Similarly, in figure 4(c), 0° < Q i < 90° for 

= —90° or 90° indicates the YZ-plane. Because 
all RCS calculations in this report are made for the 
XZ-, YZ - , and XY-planes, this notation is followed 
throughout in this report. From figure 4 we can 
conclude that M — N > 19 (which corresponds to 


the subdomain size < A 0 /10) gives reasonably stable 
results. 

Large values of M and N are required to obtain 
stable results when a flat plate is illuminated by a 
plane wave at the grazing angle because of a rapid 
variation of the phase of incident electric field along 
the plate surface. To get the quantitative estimation 
of the stability for grazing incidence, the RCS of the 
hexagonal plate shown in figure 2 is calculated as a 
function of for 6 Z = 90° and is shown in figure 5. 

The validity and accuracy of the present formu- 
lation are compared with measured results and with 
the results obtained by the Electromagnetic Surface 
Patch (ESP) code in figure 6 for parametric values 
given in the figure. For the measurement of the RCS, 
a hexagonal plate with dimensions shown in figure 2 
with a thickness equal to 0.16 cm was fabricated and 
the RCS was measured in the Langley Experimental 
Test Range facility. While using the ESP code, the 
number of modes selected was equal to 511, whereas 
the number of modes used for the present method 
was 524. Figures 3-6 show that the RCS calcula- 
tions of this paper compare reasonably well with the 
measurements and numerical results obtained from 
the ESP code. 

Figure 7 shows the RCS of the hexagonal plate 
that was calculated for grazing incidence by using the 
present formulation along with the measured data 
and results obtained using the ESP code. According 
to figure 7, the present method shows good agree- 
ment with the measured data for the incidence an- 
gles normal to the edges of the plate. However, 
for incidence angles along the corners of a plate, a 
disagreement occurs between the two results which 
may be attributed to inadequate modeling of the 
corners. 
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in figures 14 and 15, the RCS’s of this plate are 
calculated by using the present method along with 
the results obtained from the ESP code and the 
measured data. One can observe that the present 
method gives RCS estimates that are comparable 
with the measured data and with the ESP-code 
calculation . 

To compare frequency dependencies, the mono- 
static RCS’s of the plate shown in figure 13 are cal- 
culated for two angles of incidence (0* = 90° , = 90° 

and e t = 90°, & = -90°) as a function of fre- 
quency, and these data are presented in figure 16 
along with measured and calculated results from 
the ESP code. For the incidence angle 9 l = 90°, 
(fy l = 90° (the incidence on the tip), the predicted 
RCS’s are comparable with the other results. How- 
ever, for the incidence angle 6 X = 90°,^ = —90 
(the edge on incidence), the nulls and peaks pre- 
dicted by the present method are shifted in frequency 
from those observed in the measured data and in 
the results obtained from the ESP code. Note in 
these calculations that the number of modes used 
for the present method was 217 and the number of 
modes used for the ESP code was marginally smaller 
than 217. 

RCS of Equilateral Triangular Plate With 

Inverted, Equilateral Triangular J|q le 

An equilateral triangular plate with an inverted, 
equilateral triangular hole as shown in figure 17 is 
considered next. Monostatic RCS’s computed using 
the present method for various angles of incidence 
and polarization are presented in figures 18 and 19 
and compared with the results obtained from the 
ESP code. The number of modes used in both 
methods was the same. The agreement between the 
two methods was reasonably good. Because of the 
unavailability of experimental results for _this case, 
the results shown in figures 18 and 19 could not be 
compared with the measured data. 

The monostatic RCS’s of the equilateral triangu- 
lar plate with an inverted equilateral triangular hole 
as a function of frequency for a given angle of in- 
cidence are calculated by using the present method 
and are presented in figure 20 along with the results 
obtained from the ESP code. Good agreement was 
observed between the results of both methods. 


RCS of Equilateral Triangular Plate 

As a second example, an equilateral triangular 
plate with a = 5.08 cm is considered, as shown in 
figure 8. Monostatic RCS’s of the plate are calcu- 
lated as a function of the angle of incidence 9 t and 
are shown in figure 9 for parametric values given in 
this figure. From this figure, we again conclude that 
M = N > 19 gives stable results for the triangular 
plate. From these values of M and N and because 
the hexagonal and equilateral triangular plates con- 
sidered so far in this paper have dimensions around 
2Ao, one also can conclude that a subcell size of 
approximately A/ 10 is sufficient for obtaining sta- 
ble results. Also, because all the plates considered 
in this paper have maximum dimensions of approxi- 
mately 2Ao, all future calculations arc performed us- 
ing M = N = 19. However, for larger sized plates 
(larger than 2Ao), larger values of M and A are re- 
quired for convergence in accord with the A/ 10 sub- 
cell criterion. 

For parametric values given in figure 10, the RCS 
of the triangular plate is calculated by using the 
present method and is compared both with the RCS 
computed by using the ESP code and with the RCS 
measured in the Langley Experimental Test Range 
facility. The RCS of the triangular plate for grazing 
incidence is also calculated and presented in figure 11. 
From these figures we can conclude that the RCS 
predicted by the present method is in good agreement 
with both the measured data and the results obtained 
from the ESP code. Note in the results shown 
above that the number of modes selected for the 
ESP code was 312, whereas the number of modes 
P + Q for the present method was 324. Because 
of the unavailability of measured data, the results 
in figures 10(a) and 11 cannot be compared with 
measurements. 

To validate the present technique over a wide 
frequency band, monostatic RCS’s of the triangular 
plate shown in figure 8 are calculated for a fixed angle 
of incidence as a function of frequency, and these are 
presented in figure 12 along with the measured data 
and the results obtained using the ESP code. The 
results obtained by using both methods (the ESP 
code and the present method) compare very well with 
the measurements, 

RCS of Equilateral Triangular Plate With 

Concentric, Equilateral Triangular Hole 

An equilateral triangular plate with a concentric, 
equilateral triangular hole as shown in figure 13 with 
a = 5.08 cm and b = 2.54 cm is considered here as 
a third example. For the parametric values given 


RCS of Diamond-Shaped Plate 

By using the present formulation, the RCS’s of 
a diamond-shaped plate as shown in figure 21 are 
calculated and presented in figures 22 and 23 along 
with measured data and the results obtained from 
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the ESP code. Again, the agreement among the 
techniques and measurements is quite good. The 
number of modes for the two methods was again 
selected to be the same. 

Comparison of CPU Time 

From the results presented so far, the present 
method predicts the RCS of a polygonal plate with 
reasonable accuracy. The number of modes required 
to achieve this accuracy is slightly higher than the 
number of modes required by the ESP code for the 
same accuracy. This agreement is due to the fact 
that the edge conditions on the surface currents 
are not explicitly taken into account in the present 
formulation. To see the merits of the present method, 
one must compare the CPU time required to run the 
present code with the CPU time required by the ESP 
code. For a CPU time comparison, the RCS of the 
hexagonal plate shown in figure 2 is calculated by 
using the ESP code and the present method. The 
CPU time required by both methods is calculated as 
a function of the number of patch current modes and 
is plotted in figure 24. For these time comparisons, 
both codes were run on a Sun SPARCstation. 

As seen in figure 24, the present method with tri- 
angular decomposition techniques takes considerably 
less CPU time than the ESP code. This is mainly due 
to the block Toeplitz nature of the impedance matrix 
in the present method. Because the edge conditions 
on the surface currents are not explicitly taken into 
account, the present method requires more modes 
than the ESP code for reasonably accurate results. 
However, from the earlier comparisons of the results 
obtained by both methods, the difference between 
the number of modes used for the present method 
and the ESP code is marginal. (The present method 
requires 524 modes as compared with 511 modes re- 
quired by the ESP code.) The marginal increase in 
the number of modes for the present method may 
still be preferred from the CPU time consideration. 


Figure 24 also gives the CPU time in seconds as 
required by the present method using the matrix in- 
version method and the matrix triangular decompo- 
sition method. Figure 24 shows that the matrix in- 
version method takes a much longer CPU time than 
that taken by the matrix decomposition method. 

Concluding Remarks 

The method of moments has been described to 
determine the monostatic radar cross section (RCS) 
of irregularly shaped, thin, metallic flat plates in free 
space. The surface-current density on an irregularly 
shaped plate has been expressed in terms of rectangu- 
lar subdomain functions by enclosing the plate by a 
rectangle and dividing the rectangle into rectangular 
subdomains. A shape function has been used to en- 
sure zero current outside the arbitrarily shaped plate. 
The electric field integral equation (EFIE) has been 
used in conjunction with the method of moments to 
determine the surface-current density on the irregu- 
larly shaped plate. The monostatic RCS’s of (1) a 
hexagonal plate, (2) an equilateral triangular plate, 
(3) an equilateral triangular plate with a concentric, 
equilateral triangular hole and an inverted, equilat- 
eral triangular hole, and (4) a diamond-shaped plate 
have been compared with measured data and with 
the results obtained by using the Electromagnetic 
Surface Patch (ESP) code. From these comparisons, 
one can conclude that the present method calculates 
the RCS’s of these geometries with good accuracy. 
Because an irregularly shaped plate is divided into 
identical rectanglar subdomains in the present formu- 
lation, the resulting impedance matrix in the moment 
method is not only symmetrical but also of block 
Toeplitz nature. This effect is shown to result in 
a considerable time savings in filling the impedance 
matrix without sacrificing accuracy. 

NASA Langley Research Center 
Hampton, VA 23681-0001 
July 26, 1993 
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Appendix 


Description of Floppy Disk Contents 

A double-sided/high-density (DS/HD) 3j-iri. floppy disk containing several data files has been provided for 
further clarification of this report. These data files contain the experimental data and the computed results 
used in plotting various figures in the report. The computed and experimental data for a figure are stored in a 
file having the same name as the figure. For example, the data used to plot figure 3(a) may be found in a file 
named fig3a.dat. 

The floppy disk may be read using an IBM PC, PC/XT, or PC/AT with MS-DOS 2.1 or higher and a 
35-in. drive. The data on the disk are in a compressed form. To read the data, the Me figure.zip must be 
uncompressed first. This is done as shown below. Create a directory called temp , and then use the command 

pkunzip figure.zip c: temp 
This will write all data files in the temp directory. 

The reader can obtain a copy of this floppy disk by contacting the authors at the following address: 

Guidance and Control Division 
Antenna and Microwave Research Branch 
NASA Langley Research Center 
MS 490 

Hampton, VA 23681-0001 
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(a) Plate excited by plane wave. 


Y A 

m • A-directed subdomains (/?=!, 2, 3 ^ ........^P) 

O y-directed subdomains (g= 1,2, 3 } . 0) 

Approximated 

boundary of , — Actual boundary > — subdomain for 

plate — / of plate / y-directed current 



' — (m 9 /7)th subdomain for 
A-directed current 


(b) Plate divided into rectangular subdomains. 

Figure 1. Geometry of irregularly shaped, thin, metallic flat plate 




Figure 2. Thin, metallic flat plate of hexagonal shape with a = 2.074 cm. 
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(a) x- comp orient along XX- plane with //-polarized plane wave, 
x 10 2 



x!W x 


(b) y-component along XX-planc with H-polarized plane wave. 

Figure 3. Magnitude of x- and ^-components of surface- current density along XX- and yF-plane^on hexagonal 
plate (shown in fig. 2) excited by H- and JS-polarized plane waves with angle of incidence (0j - 80 , & - U ) 
for M — N = 15, 19, 21, and 23. 





(c) y-componcnt along XX-plane with ^-polarized plane wave. 


x in 2 



yw y 

(d) ./.'-component along > Y”- plane with //-polarized plane wave. 
Figure 3. Continued. 





9 /.deg 

(a) .//-polarized plane wave (a 0 = 0°) at fa — 0° and 180°. 



9 /,deg 

(b) ^-polarized plane wave (a 0 = 90°) at fa — 0° and 180°. 

Fig nre 4. Monostatic RCS of hexagonal plate (shown in fig. 2) excited by //- and ^-polarized plane waves 
(a 0 = 0 and 90 , respectively) as a function of 6 1 for / = 11.811 GHz with M = N= 15, 19, 21 and 23. 







Figure 5. Monostatic RCS of hexagonal plate (shown in fig. 2) excited by E-polarized plane wave (a 0 = 90°) 
as a function of 6, for / = 11.811 GHz with M = N = 15, 19, 21, and 23 and d l = 90° (grazing incidence). 





90 60 30 0 30 60 90 

0/.deg 

(b) ^-polarized plane wave (oq = 90°) at 4>i = 0° and 180°, 


Figure 6. Monostatic RCS of hexagonal plate (shown in fig. 2) excited by H- and E’-polarized plane waves 
(qq = 0° and 90°, respectively) as a function of 6 t for /= 11.811 GHz with M = N = 19* 






0 90 180 270 360 
<!>,■. deg 

Figure 7. Monostatic RCS of hexagonal plate (shown in fig. 2) excited by E’-polarized plane wave (a 0 = 90°) 
as a function of fa for / = 11.811 GHz with M = N = 19 and 0 t = 90° (grazing incidence). 
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(a) //-polarized plane wave (ao = 0°) at — 90° and —90 . 



(b) //-polarized plane wave (ag = 90°) at <f>i = 90° and —90°. 


Figure 9. Monostatic RCS of equilateral triangular plate (shown in fig. 8) excited by H- and E-polarizcd plane 
waves (ao = 0° and 90°, respectively) as a function of 9 l for / = 11.811 GHz with M = N = 11, 15, 19, and 



(c) H- polarized plane wave (ao = 0°) at <j>i = 0° and 180' 



6/. deg 

(d) ^-polarized plane wave (ao = 90°) at <fii = 0° and 180 
Figure 9. Concluded. 






Figure 10. Monostatic RCS of equilateral triangular plate (shown in fig. 8) excited by H- and E-polarized plane 
waves (c*o = 0° and 90°, respectively) as a function of 6 X for / = 11.811 GHz and M = N = 19. 
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Or 


Figure 11 

(c*0 = 


-10 


Present method 
ESP code 



. Monostatic RCS of equilateral triangular plate (shown in fig. 8) excited by £ , -polarized plane wave 
90°) as a function of fa for / = 11.811 GHz with M = N = 19 and 0* = 90° (grazing incidence). 
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Present method 
ESP code 
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90 60 30 0 30 60 90 

0, , deg 


(a) /?- polarized plane wave (ftQ = 90°) at — 90° and —90°. 


HH 



(b) //-polarized plane wave (ao = 0°) at fa — 90° and -90°. 

Figure 14. Monostatic RCS of equilateral triangular plate (shown in fig. 13) excited by //- and E-polarized 
plane waves (a 0 = 0° and 90°, respectively) as a function of 0; for / = 11.811 GHz with M -N - 19. 






(c) //-polarized plane wave (qq = 0°) at </>i = 0° and 180°. 



(d) /^-polarized plane wave (qq = 90°) at = 0° and 180°. 
Figure 14. Concluded. 




Figure 15. Monostatic RCS of equilateral triangular plate with concentric, equilateral triangular hole (shown 
in fig. 13) excited by /'’-polarized plane wave (ao = 90°) as a function of (p l for / = 11.811 GHz with 
M = N = 19 and = 90° (grazing incidence). 
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Frequency, GHz 

(a) 6 l = 90° , (f>i = 90° with grazing incidence at tip. 



Frequency, GHz 

(b) 6i = 90° , (j>i = —90° with grazing incidence at edge. 


Figure 16. Monostatic RCS of equilateral triangular plate with concentric, equilateral triangular hole (shown 
in fig. 13) excited by E-polarized plane waves (ao = 90°) at angles of incidence as a function of frequency 
for M= N = 19. 
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Figure 18. Monostatic RCS of equilateral triangular plate with inverted, equilateral triangular hole (shown in 
fig. 17) excited by H- and E-polarized plane waves («o = 0° and 90°, respectively) as a function of 0 t for 
/ = 11.811 GHz with M — N — 19. 
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fy.deg 

Figure 19. Monostatic RCS of equilateral triangular plate with inverted equilateral triangular hole (shown 
in fig. 17) excited by £ , -polarized plane wave (ao = 90°) as a function of (fti for / = 11.811 GHz with 
M = N = 19 and 6 Z = 90° (grazing incidence). 
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Frequency, GHz 
(a) Oi = 90°, 4>i = 90°. 



Frequency, GHz 
(b) Oi - 90°, </>i = -90°. 


Figure 20. Monostatic RCS of equilateral triangular plate with inverted, equilateral triangular hole (shown in 
fig. 17) excited by ^-polarized plane wave (ao = 90°) at angles of incidence as a function of frequency for 
M = I 




(a) H-polarized plane wave (qq — 0°) at fa = —90° and 90°. 



B /5 deg 

(b) i?-polarized plane wave (ao — 90°) at fa = —90° and 90°. 

Figure 22. Monostatic RCS of diamond-shaped plate (shown in fig. 21) excited by H- and F-polarized plane 
waves (ao = 0° and 90°, respectively) as a function of 0 t for / = 11.811 GHz with M = N = 19. 
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Figure 23. Monostatic RCS of diamond-shaped plate (shown in fig. 21) excited by E-polarized plane wave 
(ao = 90°) as a function of <f>% for / = 11.811 GHz with M = N = 19 and 0 l = 90° (grazing incidence). 



Figure 24. Comparison of CPU time required by present method using matrix inversion and decomposition 
with CPU time required by ESP code. 
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